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ABSTRACT 

The application of numerical optimization techniques to 
the design of drive shafts is demonstrated. The analysis 
investigated the effects of a small mass imbalance in con- 
junction with the rotation of a shaft with synchronous 
whirl . 
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I. 



INTRODUCTION 



Optimization and composite materials are two fields 
that have seen considerable developments in recent years. 
Advancements in these technologies have provided the design 
engineer with new dimensions, the depths of which are yet 
to be explored. 

The major material characteristics that are advanta- 
geously cultivated in composites are the high strength and 
stiffness to weight ratios. The weight reduction is so 
significant that the aircraft and aerospace industries have 
concentrated a major portion of their reserach efforts in 
this field. Investigations in improving the resistance of 
composites against environmental factors such as temperatures, 
corrosion, and wear are continuing and findings affirm the 
contention that the future of this field is bright. However, 
other industries are slow in utilizing these materials in 
basic machinery elements, mainly due to a major drawback — 
the high cost of manufacturing and processing. It is in 
this vein that the importance of an optimized design becomes 
apparent. Advancements in fabrication and processing of 
composite materials coupled with the capability to optimize 
designs may well lead composites into a more competitive 
market. 
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With the ultimate goal of optimizing the design of a 
drive shaft made of composite materials, a research program 
v/as initiated. This report is the first in this investiga- 
tive effort. 

Before the intricacies introduced by the directional 
material properties of composites are included in the 
analysis, it was deemed necessary to probe the characteris- 
tics of the design of an isotropic material. This will 
enable one to foresee the trends and the adaptive measures 
necessary when the composite material properties are con- 
sidered. The application of numerical optimization tech- 
niques to the design of an isotropic drive shaft is demon- 
strated here. The effects of a small mass imbalance on the 
bending moment and deflection of a shaft with synchronous 
whirl are investigated. Mainly, however, the objective is 
to develop an analysis procedure that may be adapted to the 
peculiarities of composites and ultimately be used in op- 
timizing a composite drive shaft design. 

Section 2 defines the scope and lim.itations of the 
analysis and optimization. 

Section 3 discusses the analysis itself, the equations 
and the modes of failure that were considered. Examples of 
the analysis are shown. 

Section 4 briefly describes the optimization program 
[1,2] and presents the results of optimizing a number of 
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shaft designs. Two variable function space diagrams are 
provided to illustrate the design space for different load 
conditions . 

Section 5 contains recommendations for future investi- 
gations . 

Appendix A contains a description of the FORTRAN program 
used in the analysis and the program text. Appendix B shows 
the derivation of equations for deflection, moment and shear 
force . 
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II. SCOPE AND LIMITATIONS 



The analysis applies to shafts made of isotropic ma- 
terials with simply supported or fixed-fixed end conditions. 
Basic shaft design formulas are used except that the de- 
flection and moment equations are derived to include the 
effects of a small mass imbalance. The shaft is ass\amed 
to rotate with synchronous whirl. 

Ta insure a continuous design field, the optimization 
is constrained to subcritical speeds only. Supercritical 
speeds introduce feasible design regions that are disjoint 
from the primary design space. In this investigation, these 
secondary areas are considered infeasible. 

Since the optimization tends to a large radius-small 
thickness design, radial stresses are neglected. Hov/ever, 
buckling of thin cylinders due to torsion and compression 
is considered. The equations used for these failure criteria 
do not include dynamic effects. A constant torque and a 
steady, uniform axial force are the loads considered in 
the formulas [3]. 

The analysis incorporates the capability to design a 
shaft that may be used in tv/o or more loading conditions. 

For example, a shaft may be designed to transmit 50 HP at 
300 RPM with an axial load of 2000 lbs. as well as to transmit 
150 HP at 3000 RPM with a 1000 lb. axial load. 
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III. ANALYSIS 



The problem is a shaft required to transmit a specified 
horsepower at a given speed. There may be an axial load, F 



are used throughout the analysis and fundamentals such as 
formulas for area, volume, moment of inertia, polar moment 
of inertia and radius of gyration are not repeated here. 
However, certain stress calculations will be emphasized to 
illustrate the flow of analysis. 

Appendix B shows the derivations of the equations for 
bending moment, deflection and shear for pinned-pinned and 
clamped-clamped end conditions. These derivations have 
taken into consideration the effects of a small mass imbalance 
in conjunction with the rotation of a shaft with synchronous 
whirl. 

The inch-pound-second system of units (IPS) is used in 
the computer program analysis. However, by changing the 
equation for torque (line 50 in the FORTRAN text. Appendix A) 
to reflect the International System of Units, the program 
may be used effectively with consistent SI units. 




(See Figure 1.) Basic equations 
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Figure 1. Applied Loads: 



Axial load (F ) , Internal pressure (P) , 
Torque (T) , Mass imbalance (e) . 




a. Pinned-pinned end condition 




b. Clamped-claraped end condition 



Figure 2. Assumed deflected shape of the shaft. 
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A. SIMPLY SUPPORTED BEAM 

The deflected shape of the beam is assumed as: 

TTX 

6=5^ sin (see Fig. 2) . The solution of the govern- 
ing eqtrdtlunS" equations is given in Appendix B. The deflection 
is maximum at x = 1/2 



max 



384 



[ 



K^e t K2 



EI-K - F(^)^ 

1 TT 7T 



] 



(3.1a) 



where: = pAw 

p = mass per unit volume 
A = cross sectional area 
w = shaft speed is radians per second 
K = pAg = weight per unit length 

e = eccentricity of mass v/ith respect to the axis 
of rotation 

a = shaft length 

E = Young's Modulus 

I = Moment of Inertia 

F = Axial Load 

Certain characteristics of this equation are of particular 
interest. If there is no rotation nor axial load, the formula 
reduces to the classical equation for maximum beam deflection. 



« 5K d 

max = 



384EI 



(3.1b) 
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If there is no axial load, the equation becomes unstable 

• 2 5, 4 

(deflection goes to infinity) v/hen El = pAco (— ) or, solving 

for the critical speed. 



(Jj = 7T 

C 



2 




(3.1c) 



which defines the first fundamental frequency of a simply 
supported shaft. 

If there is no rotation, instability arises when El = F ( 
or 



recognizable as Euler's column buckling criteria. 

With both rotation and axial load, the critical speed 
becomes 



To avoid the instability regions, the denominator of 
equation (3.1a) must be given a lower bound greater than zero 
(point A in Fig. 3) . Additionally, a maximum deflection has 
to be imposed to prevent computer overflow (point B in Fig. 3) . 
These limits have no relevance with the deflection constraint 
used in the design optimization process (point C in Fig. 3) , 
but are used only to prevent numerical ill-conditioning on 
the computer. 



F 



c 



2 

£ 



(3. Id) 




(3.1e) 



15 






Figure 3 




Limits on the deflection. The solid curve repre- 
sents the deflection values used in the analysis. 
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The bending moment corresponding to the deflection calcu- 



lated by equation (3,1a) at x = %/ ^ is 



M 



= [ I + 6 (|)^ ] + F6 + 



max 






(3. If) 



When the weight tends to cancel the deflection, the corre- 
sponding moment at the same cross section becomes 



M . = M 

min max 



4 



(3.1g) 



The shear force is maximum at x = 0. 



K,6£ K,eil K,£ 

V = — = — + — = — + — ^ 
max TT 2 2 



(3.1h) 



B . CLAMPED-CLAMPED BEAI4 

The deflected shape of the beam is assumed as; 

6 

6 = [ 1 - cos • (See Fig. 2b.) The deflection is 

maximum at x = 5-/2 • 



6 = [- 



K^e H- 



334 El - 0.746 - 96 (-) ^ F 

1 IT 



] 



(3.2a) 



As in the simply supported beam, similar observations 
regarding special loading conditions are apparent. If there 
is neither rotation nor axial load, the formula for maximum 



beam deflection results 

4 



max 384 El 



(3.2b) 
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with no axial load, instability occurs when 
384 El = 0.746 pAw^ or 



OJ 



c 



22.68 




El 

4 

pAS, 



(3.2c) 



v/hich is the first fundamental frequency for a fixed-fixed 

beam based on the assumed deflection used here. 

When there is no rotation, instability arises when 
2 

384 El = 96 or 



F 

c 



4tt^EI 



(3. 2d) 



Euler's column buckling criteria. 

With rotation and axial load, the instability occurs at 
a smaller speed than that given by equation (3.2c), 



— 22.68 




0.02533 



(3.2e) 



The bending moment at x = corresponding to the 

deflection given by equation (3.2a) is 



M 

max 



TT 



I M M 

2 7T 



+ e£ ] + f6 




(3.2f) 



A lesser bending moment occurs at the same cross section 
when the weight tends to cancel the deflection. 
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M . 
min 



M 

max 



(3.2g) 
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The maximxam shear occurs at x = 0. 



V 

max 



~ 1 ~ 



~T~ 



+ 



2 



(3.2h) 



The preceding calculations are performed in Subroutine 
BEND, given in Appendix A. 

C. STRESS ANALYSIS 

With the deflection, shear and bending moments known at 
the critical cross sections, the analysis proceeds v;ith the 
calculation of the stresses at the critical stress elements. 
At midspan, the critical stress elements are on the outer 
surface of the shaft where the axial stress (tensile or 
compressive) is maximxm (Points A and B, Fig. 4) . Since the 
shaft is in synchronous whirl, the same stress element. A, 
will generally tend to be in tension throughout the rotation, 
with the magnitude of the axial stress varying as the weight 
effect becomes alternatingly additive and subtractive. 

Element A is in the direction of the mass eccentricity from 
the center of the shaft while element B is on the outer 
surface in the opposite direction. 

Most of the following equations are found in basic design 
and applied mechanics books. References will be noted when 
necessary. 



19 . 



Axial stress due to bending moment (tension/compression) 



a (max) 

X 



M r 
_ max o 

I 



(3.3a) 



a (min) 

X 



M . r 
min o 

I 



Axial stress due to axial load (compression) 



a 

X 



£ 

A 



(3.3b) 



The axial stress due to internal pressure for thin-wall 
cylinders is; 



a 



X 




(3.3c) 



Superposition with careful regard to sign conventions 
provides the maximum and minimum tensile and compressive 
axial stresses. 

Hoop stress is primarily due to internal pressure and 
rotational inertia [4] 

Pr 

a^ = -^ (3.3d) 

a = ^ [1 + (^)(^)^] (3.3e) 

H g * 4 o 

where Y = weight per unit volume 

to = rotation in radians per second 

V = Poisson's ratio 

r. = inside radius 
1 

r = outside radius 
o 

g = acceleration due to gravity 
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The shear stress at x - JI /2 is caused only by torsion 



(shear force is zero) . 



T 



T 



Tr 
o 

J 



(3.3f) 



where = shear stress due to torsion 

T = torque 

J = polar moment of inertia 

At X = 0, the critical element is where the shear stress 
due to the shear force along the cross section is additive 
with respect to the shear stress due to torsion (element C, 
Fig. 3) . The shear stress due to for thin wall cylinders 

[5] is 



2 V 



"v = 



max 



(3.3g) 



Superposition of equations (3.3f) and (3.3g) gives the total 
shear stress at the critical stress element. 

The axial stress at element C due to bending moment is 
zero. The axial stress at this point is, therefore, caused 
only by the axial load and the internal pressure, equations 
(3.3b) and (3.3c) . 

The hoop stress is calculated by equations (3.3d) and 
(3.3e) . 
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stress element A 




a) At midspan 




b) At X = 0 




c) Typical stress element. Radial 
stresses are neglected. 

Figure 4. Critical Stress ' Elements 
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D. MODES OF FAILURE 



Eight modes of failure are considered in this investiga- 
tion. The first two evaluate the static strengths at the 
critical stress elements using Von Mises criteria [5] . At 
midspan, 



Failure is assumed if the yield strength of the material is 
equal to or less than these Von Mises stresses multiplied by 
a safety factor. 

A limit on the deflection is imposed. A common practice 
is to specify a maximum deflection in inch per foot of shaft 
length. For example, as used in this report, 

6 = 0.005 inch/foot length 

max 

Since the optimization may lead to a thin cylinder design, 
buckling due to torsion and compression are considered. The 
critical torque and critical compressive stress for buckling 
of thin cylinders [3] are; 




(3.4a) 



At X = 0 




(3.4b) 



T 




(3.4c) 
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a 



Et 



(3.4d) 



r^3(l-v2)' 

Experimental data indicate that actual failure of the cylinders 
usually occurs below fifty percent of the values calculated 
from these equations. To accommodate this discrepancy and 
the fact that the dynamic effects are not considered in the 
derivation of these equations, comparatively higher safety 
factors must be used for these failure modes. 

The shaft is also checked for column buckling using Johnson's 
or Euler's equations [5] as applicable. The actual slenderness 
ratio is compared with the slenderness ratio where transition 
occurs to determine which formulas apply. The calculated 
critical load, using the appropriate equations for the end 
condition, is compared to the actual load, again incorporating 
a safety factor. Equations (3. Id) or (3. 2d) are used if the 
slenderness ratio falls in the Euler range. In the Johnson 
region, the following equations [5] apply. 



F 

c 

A 



S -b 
c 



(-) ^ 



(3.4e) 



where S = yield strength in compression 



(^) 



slenderness ratio 
S 



'■2-n^ nE 



b 

E = Young's Modulus 
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n = 1 for simple supports 
4 for fixed-fixed ends 

The effects of the fluctuating axial stresses are evaluated 
for fatigue failure using the Goodman diagram [5]. Von Mises 
stresses for the mean and alternating stresses are calculated 
as follows: 

CT + a . 

a = (3.4f) 

m 2 



a 

a 



a 

max - 



2 



a . 

min 



(3.4g) 



a 



m 




H m 



+ a 



m 



+ 3 




(3.4h) 



a ' = a (3.4i) 

a a 

The Goodman diagram is drawn (Fig. 5) and the safety factor 
is determined analytically. 




Figure 5. Goodman Diagram 
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SF 



X 



Se 



(3.4j) 



m 



m 



- + 1 ^] 
, Su-* 



m 



where Se = fully corrected endurance limit 
Su = ultimate strength 

This safety factor is required to be greater than a specified 
value. 

Lastly, failure is assumed if a specified percentage of 
the critical speed as calculated by equations (3.1e) or (3.2e) 
is less than the shaft speed. 

E. ANALYSIS EXAMPLES 

The examples are chosen to illustrate several aspects of 
the shaft design problem. The first is a design which vio- 
lates almost all of the constraints while the second example 
is one within the feasible region of the design field. 

The material used in all the examples in this report is 
UN G10150 HR steel, the properties of which are listed in 
Table I [5] . 
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TABLE I 



MATERIAL PROPERTIES 
UNS G10150 HR STEEL 



Young ' s Modulus 
Shear Modulus 
Poisson's Ratio 
Weight Per Unit Volume 
Yield Strength 
Ultimate Strength 



30. OE 06 psi 
11. 5E 06 psi 
0.292 ^ 

0.282 Ibs/in-^ 
27. OE 03 psi 
50. OE 03 psi 
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1 . Analysis Example 1 



A simply supported shaft with 0.1 inch eccentricity 
is to transmit 150 HP at 3000 RPM (3150 in-lb torque) , 
with an axial load of 2000 lbs. For a shaft length of 

120 inches, inside radius of 1.0 inch and thickness of 
1.0 inch, the design violates the static and dynamic 
strength, deflection, buckling due to compression, and 
the speed constraints. The analysis of this design is given 
below where an asterisk (*) denotes violated constraints on 
the design. 



Volume 


0.11310E 


04 


in^ 


Weight 


0. 31893E 


03 


lbs 


Moment of Inertia 


0.11781E 


02 


in^ 


Polar Moment of Inertia 


0.23562E 


02 




Critical Speed 


0.14767E 


04 


RPM 


Critical A.xial Load 


0.18764E 


06 


lbs 


Critical Buckling Torque 


0.47469E 


08 


in lb 


Critical Buckling Stress 


0.905493 


07 


psi 


Deflection 


0.24000E 


02 


in 


Bending Moment 


0. 36420E 


08 


in lb 


Shear Force 


0.62722E 


06 


lb 


Axial Stress at A 


0.95779E 


07 


psi 


Axial Stress at B 


-0.95783E 


07 


psi 


Hoop Stress 


0.25013E 


03 


psi 


Torsional Shear Stress 


0. 26738E 


03 


psi 


Total Shear Stress 


0.13337E 


06 


psi 


Mean Axial Stress 


0. 95771E 


07 


psi 


Alternating Axial Stress 


0.81250E 


03 


psi 


Von Mises at Midspan 


0.95785E. 


07 


psi 


Von Mises at x = 0 


0.23100E 


06 


psi 
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The calculated deflection is the maximum value 



imposed to avoid numerical ill-conditioning (point B, Fig. 3) . 
Since the deflection affects most of the relevant parameters 
in the analysis, the calculated stresses are unreliable and 
not valid. These values, however, tend to direct the design 
towards the feasible region. 

2 . Analysis Example 2 

This shaft is similar to the first example in design 
conditions and loading. The only difference is that the 
dimensions are such that no constraints are violated, but 
two constraints are active or critical. The shaft is 
120 inches long with an inside radius of 6.39 inches and 
thickness of 0.0369 inch. This, in fact, is the optimized 
design for the given loading as calculated by COPES-CONMIN 
[1,2]. The results of the analysis follows, where a double 
asterisk (**) denotes values at their design limit. 



Volume 


0.17866E 


03 


. 3 
in 


Weight 


0.50382E 


02 


lbs 


Moment of Inertia 


0. 30629E 


02 




Polar Moment of Inertia 


0.61258E 


02 




*Critical Speed 


0. 60061E 


04 


RPM 


Critical Axial Load 


0.39557E 


05 


lbs 


*Critical Buckling Torque 


0. 31490E 


05 


in lb 


Critical Buckling Stress 


0.10400E 


06 


psi 


Deflection 


0.43813E- 


•01 


in 


Bending Moment 


0.17026E 


05 


in lb 


Shear Force 


0.84895E 


03 


lbs 


Axial Stress at A 


0.22326E 


04 


psi 


Axial Stress at B 


-0.49193E 


04 


psi 


Hoop Stress 


0.29778E 


04 


psi 


Torsional Shear Stress 


0. 33079E 


03 


psi 


Total Shear Stress 


0.14712E 


04 


psi 


Mean Axial Stress 


0.20739E 


04 


psi 


Alternating Axial Stress 


0.15872E 


03 


psi 


Von Mises at Midspan 


0. 69313E 


04 


psi 


Von Mises at x = 0 


0. 39193E 


04 


psi 
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In each of the examples, the analysis proceeds as 
follows. All shaft characteristics like cross sectional 
area, volume, weight, moment of inertia and polar moment 
of inertia are first calculated. Parameters that depend 
only on geometry such as critical torque and compressive 
stress for buckling of thin cylinders, and critical speed 
are then determined before Subroutine BEND, where the de- 
flection, bending moments and shear force at the critical 
cross sections are computed, is called. 

At midspan, the absolute values of the axial 
stresses at elements A and B (see Fig. 4a) are compared 
to determine which element is critical. In both examples, 
element B is considered critical for strength evaluation. 

At x = 0, the critical element, C, is where the torsional 
shear stress and the shear stress due to the shear force 
along the cross section are additive. The Von Mises 
stresses for each of the critical elements are then calcu- 
lated using equations (3.4a) and (3.4b). The rest of the 
analysis is a straight forward application of the equations 
for the critical parameters and the constraints. 
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IV. OPTIMIZATION, RESULTS AND CONCLUSIONS 



A. OPTIMIZATION AND RESULTS 

There are a number of optimization programs available, 
each one using different techniques in locating the desired 
optimum design. COPES-CONMIN [1,2] is a versatile program 
that may be used for sensitivity analysis and two variable 
function space study as well as an optimization tool. 

Provided with a user supplied analysis program (Subroutine 
ANALIZ) where the objective function , constraints and other 
relevant parameters are calculated, it determines a usable 
and feasible sector from which a search direction is chosen. 
This choice of search direction is made by using such informa- 
tion as the gradients of the objective function, gradients of 
active and violated constraints and the degree of push-off 
these constraints violations generate. This iterative process 
of minimizing/maximizing the objective function by changing 
the design variables is terminated when no further improvement 
can be made. 

In this exercise, the objective function to be minimized 
is the volume, with the radius and thickness as design 
variables. Several design examples are shown to demonstrate 
the effects of the different parameters on the optimized 
design and on the design space itself. Seven examples are 
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presented corresponding to the loading conditions as listed 
in Table II. Two variable function space diagrams are 
generated to illustrate these effects. The numbered curves 
in the graphs represent the corresponding constraints shown 
in Table III. The shaded side of a curve represents the un- 
feasible area of the design space with respect to the cor- 
responding constraint. The dashed curves represent contours 
of the objective function, the numbers indicating the 
volume in cubic inches. 



TABLE II 

OPTIMIZATION EXAMPLES 



No. 


Loads 


Horsepower 


RPM 


Pressure 

(psi) 


Axial Load 
(lbs) 


1 


150 


300 


0 


0 


2 


150 


300 


1000 


0 


3 


150 


300 


0 


2000 


4 


150 


3000 


0 


0 


5 


150 


3000 


1000 


0 


6 


150 


3000 


0 


2000 


7* 


150 


300 


0 


0 




150 


3000 


0 


2000 



*Example No. 7 is a design for two load conditions. 
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TABLE III 



DESIGN CONSTRAINTS 



Nvunber 



Constraint 



1 

2 

3 

4 

5 

6 

7 

8 



Static Strength at Midspan 
Static Strength at x = 0 
Deflection 

Buckling due to Torque 
Buckling due to Compression 
Column Buckling 
Dynamic Strength (Fatigue) 
Speed 
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EX^^PLE 1. : INITIAL OSSISM 

SHAFT DIMENSIDNS : 

RI = O.IOOOOE 01 TH = O.IOOOOE 01 SL = 0.12000E 
DESIGN CONDITIONS : 

ECCENTRIGITT = 0.100005 00 
END CONDITIONS : P INNED-PI NNED 
LOADS : 

HORSEPOWER RPM PRESSURE AXIAL LOAD 

0.15000E 03 0.30000E 03 0.0 0.0 



ANALYSIS RESULTS 



TORQUE! U 


0.31500E 


05 


VOLU'IE 


0. 1131DE 


04 


WEIGHT 


0.31893E 


03 


MOMENT OF INERTIA 


0.11731E 


02 


POLAR MOMENT OF INERTIA 


0.23562E 


02 


CRITICAL SPEED 


0.14828E 


04 


CRITICAL AXIAL LOAD 


0.18754E 


06 


CRITICAL BUC<LING torqjc 


0.47<»59E 


08 


CRITICAL BUC<LING STRESS 


0.905^95 


07 


DEFLECTION 


0.26584E- 


-01 


BENDING MOMENT 


0.50477E 


04 


SHEAR FORCE 


0.2071iVE 


03 


AXIAL STRESS AT A 


0.856926 


03 


AXIAL STRESS AT B 


-0.856925 


03 


HOOP STRESS 


0.25013E 


01 


TORSIONAL SHEAR STRESS 


0.26738E 


04 


TOTAL SHEAR STRESS 


0.27178E 


04 


MEAN AXIAL STRESS 


0.44760E 


02 


ALTERNATING AXIAL STRESS 


0.81216E 


03 


VON MISES AT MIDSPAN 


0.47095E 


04 


VON MISES AT X=0. 


0.47073E 


04 
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EXAMPLE 1. : OPTIMUM OESIG'< 

SHAFT DIMENSIONS : 

RI = 0.25847= 01 TH = 0.U124E 00 SL 
DESIGN CONDITIONS : 

ECCENTRICIT7 = O.IOOOOE 00 
END CONDITIONS : P I NNEO-PI NNED 
LOADS : 

HORSEPOWER RPM PRESSURE 

0.15000E 03 0.30000E 03 0.0 



= 0.12000E 



AXIAL LOAD 
0 . 0 



ANALYSIS RESULTS 



TORQ'JEd) 


0.31500E 


05 


VOLUME 


0.221455 


03 


WEIGHT 


0.62451E 


02 


MOMENT OF INERTIA 


0.64353E 


01 


POLAR MOMENT OF INERTIA 


0.123715 


02 


CRITICAL SPEED 


0.247575 


04 


CRITICAL AXIAL LOAD 


0.45137E 


05 


CRITICAL BUC<LING TORQUE 


0.315005 


06 


CRITICAL BUCKLING STRESS 


0. 74729 E 


06 


DEFLECTION 


0.92753E- 


-02 


BENDING MOMENT 


0.95479E 


03 


SHEAR FORCE 


0.39681E 


02 


AXIAL STRESS AT A 


0.39998E 


03 


AXIAL STRESS AT B 


-0.39998E 


03 


HOOP STRESS 


0.51655E 


01 


TORSIONAL SHEAR STRESS 


0.65980E 


04 


TOTAL SHEAR STRESS 


0.664105 


04 


MEAN AXIAL STRESS 


0.755075 


01 


ALTERNATING AXIAL STRESS 


0.39243E 


03 


VON MISES AT MIDSPAN 


0.11435E 


05 


VON MI SES AT X=0. 


0.11503E 


05 
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inside; RNDIU5 INCHtS C 



EXVMPLE 2 



IMITIAL OESIGVJ 



SHAFT DIMENSI3MS : 

RI = 0.10000= 01 TH = O.IOOOOE 01 SL = 0.12000E 
DESIGN CONOITI0NS : 

ECCENTRICITY = O.IOOOOE 00 
END CONDITIONS : PI NNED-PI NN5D 
LOADS : 

HORSEPOWER RPM PRESSURE AXIAL LOAD 

0.15000E 03 0.30000E 03 O.IOOOOE 0^ 0.0 



ANALYST S RESULTS 



TORQJE(l) 


0.31500E 


05 


VOLUME 


0. 11310E 


04 


WEIGHT 


0.31893E 


03 


MOMENT OF INERTIA 


0.11781E 


02 


POLAR MOMENT OF INERTIA 


0.23562E 


02 


CRITICAL SPEED 


0. 14828E 


04 


CRITICAL AXIAL LOAD 


0. 1876AE 


06 


CRITICAL BUC<LING TORQUE 


0.47A69E 


08 


CRITICAL BUC<LING STRESS 


0.90549E 


07 


DEFLECTION 


0.26584E- 


•01 


BENDING MOMENT 


0.50477E 


04 


SHEAR FORCE 


0.2071^E 


03 


AXIAL STRESS AT A 


0.18569E 


04 


AXIAL STRESS AT B 


0.14308E 


03 


HOOP STRESS 


0.66917E 


03 


TORSIONAL SHEAR STRESS 


0.26738E 


04 


TOTAL SHEAR STRESS 


0.27178E 


04 


MEAN AXIAL STRESS 


0.104^8E 


04 


ALTERNATING AXIAL STRESS 


0.81216E 


03 


VON MI SES AT MIDSPAN 


0.49093E 


04 


VON MI SES AT X = 0. 


0.475A&E 


04 
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EXXMPLE 2. : DPTI'^UM OESIGvj 

SHAFT DIMENSIONS : 

RI = 0.18625E 01 TH = 0.2S560E 00 5L = 0.12000E 
DESIGN CONDITIONS : 

ECCENTRICITY * O.IOOOOE 00 
END CONDITIONS : PI NNED-PI NNED 
LOADS : 

HORSEPOWER RPM PRESSURE AXIAL LOAD 

0.15000E 03 0.30000E 03 O.IOOOOE 0^ 0.0 



ANALYSIS RESULTS 



TORQJEd ) 


0.31500E 


05 


VOLUME 


0.39958E 


03 


WEIGHT 


0.11268E 


03 


MOMENT OF INERTIA 


0.66578E 


01 


POLAR MOMENT OF INERTIA 


0.13316E 


02 


CRITICAL SPE=0 


0.1875^E 


04 


CRITICAL AXIAL LOAD 


0.751^4E 


05 


CRITICAL BUC<LING TORQUE 


0.23553E 


07 


CRITICAL BUC<LING STRESS 


0.22602E 


07 


DEFLECTION 


0.16358E- 


•01 


BENDING MOMENT 


0.17476E 


04 


SHEAR FORCE 


0.72247E 


02 


AXIAL STRESS AT A 


0.75709E 


04 


AXIAL ■ STRESS AT B 


0.64537E 


04 


HOOP STRESS 


0.65487E 


04 


TORSIONAL SHEAR STRESS 


0.50344E 


04 


TOTAL SHEAR STRESS 


0.50773S 


04 


MEAN AXIAL STRESS 


0.70306E 


04 


ALTERNATING AXIAL STRESS 


0.54027E 


03 


VON MISES AT MIOSPAN 


0.11254E 


05 


VON MISES AT X=0. 


0.10965E 


05 
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inside; NBOIUS INCHETS 



ex^'^PLE 3 



nlTIAL 0ESI3N 



SHAFT DIMENSIONS : 

RI = 0.100005 01 TH = O.IOOOOE 01 SL = 0.12000E 
DESIGN CONDITIONS : 

ECCENTRICITY = O.IOOOOE 00 
END CONDITIONS : PI NNED-PI NNED 
LOADS : 

HORSEPOWER RPM PRESSURE AXIAL LOAD 

0.15000E 03 0.30000E 03 0.0 0.20000E 04 



ANALYSIS RESULTS 



TORQUE m 


0.31500E 


05 


VOLUME 


0.11310E 


04 


WEIGHT 


0.31893E 


03 


MOMENT OF INERTIA 


0.11731E 


02 


POLAR MOMENT OF INERTIA 


0.23562= 


02 


CRITICAL SPE50 


0. 14767E 


04 


CRITICAL AXIAL LOAD 


0.18754E 


06 


CRITICAL BUCXLING TORQUE 


0.47469E 


08 


CRITICAL BUC<LING STRESS 


0.90549E 


07 


DEFLECTION 


0.26814E- 


•01 


BENDING MOMEviT 


0.54144E 


04 


SHEAR FORCE 


0.20720E 


03 


AXIAL ■ STRESS AT A 


0.70697E 


03 


AXIAL STRESS ^T B 


-0. 11314E 


04 


HOOP STRESS 


0.25013E 


01 


TORSIONAL SHEAR STRESS 


0.26739E 


04 


TOTAL SHEAR STRESS 


0.27178E 


04 


MEAN AXIAL STRESS 


-0.10519E 


03 


ALTERNATING AXIAL STRESS 


0.81216E 


03 


VON MISES AT MIDSPAN 


0.47677E 


04 


VON MISES AT X=0. 


0.470735 


04 
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EXAMPLE 3. : OPTIMUM DESIGM 

SHAFT DIMENSIONS : 

RI = 0.259075 01 TH = 0.11119E 00 SL = 0.12000E 
DESIGN CONDITIONS : 

ECCENTRICITY = O.IOOOOE 00 
END CONDITIONS : PI NNEO-PI NN50 
LOADS : 

HORSEPOWER RPM PRESSURE AXIAL LOAD 

O.15DO0E 03 0.30000E 03 0.0 0.20000E OA- 



ANALYSIS RESULTS 



TOROUE(l) 


0.31500E 


05 


VOLUME 


0.22186E 


03 


WEIGHT 


0.6256!»E 


02 


MOMENT OF INERTIA 


0.64-766E 


01 


POLAR MOMENT OF INERTIA 


0.12953E 


02 


CRITICAL SPEED 


0.2463SE 


04 


CRITICAL AXIAL LOAD 


0.45241E 


05 


CRITICAL BUCCLING TORQUE 


0.315005 


06 


CRITICAL BUC<LING STRESS 


0.74527E 


06 


DEFLECTION 


0.93752E- 


■02 


BENDING MOMENT 


0.981715 


03 


SHEAR FORCE 


0.39758E 


02 


AXIAL STRESS AT A 


-0.67221E 


03 


AXIAL STRESS AT B 


-0. 14913E 


04 


HOOP STRESS 


0.51838E 


01 


torsional SHEAR STRESS 


0.65705E 


04 


TOTAL SHEAR STRESS 


0.66136E 


04 


MEAN AXIAL STRESS 


-0.10637E 


04 


ALTERNATING AXIAL STRESS 


0.39151E 


03 


VON MISES AT MIDSPAN 


0.114785 


05 


VON MISES AT X=0. 


0.11455E 


05 
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IN^lDe RRDIUS INCHETS 



SX\MPL6 ^ 



INITIAL DESIGN 



SHAFT DIM6NSI3NS : 

RI = 0.10000= 01 TH = O.IOOOOE 01 SL = 0.12000E 
DESIGN CONDITIONS : 
eCCENTRICITT = 0.100005 00 
END CONDITIONS : PI NNED-P I NNED 
LOADS : 



HORSEPOWER RPN 


PRESSURE 


AXIAL LOAD 


0.15000E 03 0.30000E 04 


o 

• 

o 


o 

• 

o 


ANALYSTS RESULTS 




TORQUE (1) 


0.31500E 


04 


VOLUNE 


0.11310E 


04 


WEIGHT 


0.31893E 


03 


MOHENT OF INERTIA 


0. 11781E 


02 


POLAR M0H5NT OF INERTIA 


0.23562E 


02 


CRITICAL SPEED 


0.14828E 


04 


CRITICAL axial LOAD 


0.18764E 


06 


CRITICAL 3UC<lING TORQUE 


0.47459E 


08 


CRITICAL 3UC<LING STRESS 


0-90549E 


07 


DEFLECTION 


0.24000E 


02 


3EN0ING MOMENT 


0.23301E 


08 


SHEAR FORCE 


0.62722E 


06 


AXIAL STRESS AT A 


0.40406E 


07 


AXIAL STRESS AT 3 


-0.40406E 


07 


HOOP STRESS 


0.25013E 


03 


TORSIONAL SHEAR STRESS 


0.26733E 


03 


TOTAL SHEAR STRESS 


0.13337E 


06 


MEAN AXIAL STRESS 


0.40397E 


07 


ALTERNATING AXIAL STRESS 


0.81200E 


03 


VON MISES AT MIDSPAN 


0.40404E 


07 


VON MISES AT X=0. 


0.23100E 


06 



43 



EXAMPLE 4 



3PTIMUM OESIGM 



SHAFT OIMENSIDNS : 

RI = 0.63798E 01 TH = 0.36963E-01 SL = 0. 12000E 
DESIGN CONDITIONS : 

ECCENTRICITY = O.IOOOOS 00 
END CONDITIONS : PI NNEO-PI NNEO 
LOADS : 

HORSEPOWER RPM PRESSURE AXIAL LOAD 

0.15000E 03 0.30000E 04 0.0 0.0 



ANALYSIS RESULTS 



TORQUEU) 


0.31500E 


04 


VOLUME 


0.17831E 


03 


WEIGHT 


0.50283E 


02 


MOMENT OF INERTIA 


0.30415E 


02 


POLAR MOMENT OF INERTIA 


0.60830E 


02 


CRITICAL SPEED 


0.60005E 


04 


CRITICAL axial LOAD 


0.39476E 


05 


CRITICAL BUC<LING TORQUE 


0.31494E 


05 


CRITICAL BUC<LING STRESS 


0.10432E 


06 


DEFLECTION 


0.43922E- 


-01 


BENDING MOMENT 


0.76214E 


04 


SHEAR FORCE 


0.84773E 


03 


AXIAL STRESS AT A 


0. 16079E 


04 


AXIAL STRESS AT B 


-0.16079E 


04 


HOOP STRESS 


0.29628E 


04 


TORSIONAL SHEAR STRESS 


0.33228E 


03 


TOTAL SHEAR STRESS 


0.14733E 


04 


MEAN AXIAL STRESS 


0.14488E 


04 


ALTERNATING AXIAL STRESS 


0.15913E 


03 


VON MISES AT MIDSPAN 


0.26327E 


04 


VON MISES AT X=0. 


0.391035 


04 
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INStOE: RHDIU5 INCHES ( 



EX^'1PLE 5. : INITUL DESIGN 

SHAFT DIMENSIDNS : 

RI = 0.100005 01 TH = O.IOOOOE 01 SL = 0. 12000E 
DESIGN CONDITIONS : 

ECCENTRICITY = 0.10000= 00 



END CONDITIONS : PINNED-PI 
LOADS : 

HORSEPOWER RPN 

0.15000E 03 0.30000E OA- 

ANALYSIS RESU 

TORQUE (U 

VOLUME 

WEIGHT 

MOMENT OF INERTIA 
POLAR MOMENT OF INERTIA 
CRITICAL SPEED 
CRITICAL AXIAL LOAD 
CRITICAL BUC<LING TORQUE 
CRITICAL 8UC<LIiNG STRESS 
DEFLECTION 
BENDING MOMENT 
SHEAR FORCE 
AXIAL STRESS AT A 
AXIAL STRESS AT B 
HOOP STRESS 
TORSIONAL SHEAR STRESS 
TOTAL SHEAR STRESS 
MEAN AXIAL STRESS 
ALTERNATING AXIAL STRESS 
VON MISES AT MI DSPAN 
VON MISES AT X=0. 



NNED 



PRESSURE 


AXIAL LOAD 


O.IOOOOE 04- 


o 

• 

o 


TS 


0.31500E 


04 


0. 11310E 


04 


0.31893E 


03 


0.11731E 


02 


0.23552E 


02 


0.14823E 


04 


0.1876^E 


06 


0.47459E 


08 


0.905A9E 


07 


0.24000E 


02 


0.23801E 


08 


0.62722E 


06 


0.40416E 


07 


-0.40396E 


07 


0. 91680 = 


03 


0.26738E 


03 


0.13337E 


06 


0.40407E 


07 


0.81200E 


03 


0.40411E 


07 


0.231005 


06 



46 



EX\>^PLE 5 



OPTIMUM 0E5IGM 



SHAFT DIMENSIONS : 

RI = 0.6084AE 01 TH = 0.6L266E 00 SL = 0. 12000E 
DESIGN CONDITIONS : 

ECCENTRICITY = O.IOOOOE 00 
END CONDITIOMS : P INNEO-PI NNEO 
LOADS : 

HORSEPOWER RPM PRESSURE AXIAL LOAD 

0.15000E 03 0.30000E 04 O.IOOOOE 04 0,0 



ANALYSIS results 



TOROUEdJ 


0.31500E 


04 


VOLJNE 


0.29521E 


04 


WEIGHT 


0.83250E 


03 


MOMENT OF INERTIA 


0.50353E 


03 


POLAR MOMENT OF INERTIA 


0. 10071E 


04 


CRITICAL SPEED 


0.60003E 


04 


CRITICAL AXIAL LOAD 


0.65357E 


06 


CRITICAL BUCKLING TOROJE 


0.34401E 


08 


CRITICAL BUCKLING STRESS 


0.16557E 


07 


DEFLECTION 


0.43926E- 


•01 


BENDING MOMENT 


0.12619E 


06 


SHEAR FORCE 


0.14035E 


05 


AXIAL STRESS AT A 


0.11609E 


05 


AXIAL STRESS AT B 


0.82525E 


04 


HOOP STRESS 


0.12589E 


05 


TORSIONAL SHEAR STRESS 


0.20948E 


02 


TOTAL SHEAR STRESS 


0.11620E 


04 


MEAN AXIAL STRESS 


0.11443E 


05 


ALTERNATING AXIAL STRESS 


0.16609E 


03 


VON MIS5S AT MIDSPAN 


0.12129E 


05 


VON MI SES AT X=0. 


0.12749E 


05 
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INSIPE; RFID1U5 INCHCS 



example 6. : INITIAL DESIGN 

SHAFT DIMENSIONS : 

RI = 0.10000= 01 TH = O.IOOOOE 01 SL = 0.12000E 
DESIGN CONDITIONS : 

ECCENTRICITY = O.IOOOOE 00 



END CONDITIONS : PINNED-PI 
LOADS : 

HORSEPOWER RPM 

0.15000E 03 0.30000E 04 

ANALYSIS RES 

TORQjEd ) 

VOLUME 
WEIGHT 

MOMENT OF INERTIA 
POLAR MOMENT OF INERTIA 
CRITICAL SPEED 
CRITICAL AXIAL LOAD 
CRITICAL 8UC<LING TORQUE 
CRITICAL BUCKLING STRESS 
DEFLECTION 
BENDING MOMENT 
SHEAR FORCE 
AXIAL STRESS AT A 
AXIAL STRESS AT B 
HOOP STRESS 
TORSIONAL SHEAR STRESS 
TOTAL SHEAR STRESS 
MEAN AXIAL STRESS 
alternating axial STRESS 
VON MISES AT MIOSPAN 
VON MI SES AT X=0. 



NNED 



PRESSURE 


AXIAL LOAD 


o 

• 

o 


0.20000E 04 


TS 


0.31500E 


04 


O.11310E 


04 


0.31893E 


03 


0. 11731E 


02 


0.235S2E 


02 


0. 14757E 


04 


0. 18764E 


06 


0.47469E 


08 


0.90549E 


07 


0.2400DE 


02 


0.56420E 


08 


0.62722E 


06 


0.95779E 


07 


-0.95783E 


07 


0.25013E 


03 


0.26733E 


03 


0.13337E 


06 


0.95771E 


07 


0.81250E 


03 


0.95785E 


07 


0.23100E 


06 
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EXAMPLE 6. : DPTI^UM DESIGVJ 

SHAFT DIMENSIONS : 

RI = 0.63959E 01 TH = 0.369^2E-0l SL = 0.12000E 
DESIGN CQNDiriDNS : 

ECCENTRICITY = 0.19000E 00 
END CONDITIONS : P! NNEO-PI NNED 
LOADS : 

HORSEPOWER RPM PRESSURE AXIAL LOAD 

0.15000E 03 0.30000E OA- 0.0 0.20000E 04 



ANALYSIS RESULTS 



TORQJE(l) 


0.31500E 


04 


VOLUME 


0.178S6E 


03 


WEIGHT 


0.50382E 


02 


MOMENT OF INERTIA 


0.306295 


02 


POLAR MOMENT OF INERTIA 


0.61253E 


02 


CRITICAL SPEED 


0.60061E 


04 


CRITICAL AXIAL LOAD 


0.39557E 


05 


CRITICAL BUC<L!NG TORQUE 


0.31490E 


05 


CRITICAL BUCXLING STRESS 


0.10400E 


06 


DEFLECTION 


0.43813E- 


•01 


BENDING MOMENT 


0.17025E 


05 


SHEAR FORCE 


0.848955 


03 


AXIAL STRESS AT A 


0.22326E 


04 


AXIAL ■ STRESS 0 


-0.49193E 


04 


HOOP STRESS 


0.29778E 


04 


TORSIONAL SHEAR STRESS 


0.33079E 


03 


TOTAL SHEAR STRESS 


0.14712E 


04 


MEAN AXIAL STRESS 


0.20739E 


04 


ALTERNATING AXIAL STRESS 


0. 15372E 


03 


VON MISeS AT MIDSPAN 


0.69313E 


04 


VON MISES AT X=0. 


0.39193E 


04 
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INSIDe RHDIUS INCHE!5 ( 



EXV'IPLE 7 

SHAFT DIMENSIONS : 



INITIAL DESIGN 



RI = O.IOOOOE 01 TH = 0 


.lOOOOE 01 SL 


= 0. 12000E 


DESIGN CONDITIONS : 








ECCENTRIGITy = 0.100005 


30 






END CONDITIONS : PINNED-PI NNED 






LOADS : 








HORSEPOWER RPM 


PRESSURE 


AXIAL LOAD 


0.15000E 03 0.30000E 03 


0.0 




o 

• 

o 


0.15000E 03 0.30000E 0^ 


0.0 




0.20000E 04 


ANALYSIS RESULTS 






TORQUE <U 




0.31500E 


05 


T0R0JE(2 » 




0.31500E 


04 


VOLUME 




O.11310E 


04 


WEIGHT 




0.31893E 


03 


MOMENT OF INERTIA 




0.11731E 


02 


POLAR MOMENT OF INERTIA 




0.23562E 


02 


CRITICAL SPEED 




0.14767E 


04 


CRITICAL AXIAL LOAD 




0.1876^E 


06 


CRITICAL BUC<LING TORQUE 




0.^7469E 


08 


CRITICAL BUCKLING STRESS 




0.90549E 


07 


DEFLECTION 




0.24000E 


02 


BENDING MOMENT 




0.56^20E 


08 


SHEAR FORCE 




0.62722= 


06 


AXIAL STRESS AT A 




0.95779E 


07 


AXIAL ■ STRESS AT B 




-0.95783E 


07 


HOOP STRESS 




0.25013E 


03 


TORSIONAL SHEAR STRESS 




0.26733E 


03 


TOTAL : SHEAR STRESS 




0.13337E 


06 


MEAN AXIAL STRESS 




0.95771E 


07 


ALTERNATING AXIAL STRESS 




0.81250E 


03 


VON MI SES AT MIDSPAN 




0.95785E 


07 


VON MISES AT X=0. 




0.23100E 


06 
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EXX^»LE 7. : DPTI'IUM 0ESI3v| 

SHAFT 0IMENSI3NJS : 



RI = 0.63551E 01 TH = 0 


.93009E-01 SL 


= 0.12000E 


DESIGN CDNOITIDNS : 








ECCENTRICITY = 0.10000= 


00 






END CONDITIONS : P INNED-PI NNED 






LOADS ; 








HORSEPOWER RPM 


PRESSURE 


AXIAL LOAD 


0.15000E 03 0.30000E 03 


o 

• 

o 




o 

• 

o 


0.15000= 03 0.30000E 04 


0.0 




0.20009E 04 


\N4LYSIS results 






TORQUE(l) 




0.31500E 


05 


TORQUE (2) 




0.31500E 


04 


VOLUME 




0.44893E 


03 


WEIGHT 




0.12660E 


03 


MOMENT OF INERTIA 




0.76660E 


02 


POLAR MOMENT OF INERTIA 




0.15332E 


03 


CRITICAL SPEED 




0.60000E 


04 


critical axial load 




0.99390E 


05 


CRITICAL BUCXLING TORQUE 




0.31571E 


06 


critical BUCXLING STRESS 




0.2S122E 


06 


DEFLECTION 




0.43932E- 


•01 


BENDING MOMENT 




0.42893F 


05 


SHEAR FORCE 




0.21344E 


04 


AXIAL STRESS AT A 




0.30733E 


04 


AXIAL STRESS AT B 




-0.41425E 


04 


HOOP STRESS 




0.29828F 


04 


TORSIONAL SHEAR STRESS 




0.13248E 


03 


TOTAL SHEAR STRESS 




0.12736E 


04 


MEAN AXIAL STRESS 




0.29136E 


04 


ALTERNATING AXIAL STRESS 




0.15973E 


03 


VON MI SES AT MIDSPAN 




0.62021E 


04 


VON MI SES AT X=0. 




0.37099E 


04 
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INSiOC RRDIUS INCHC5 




B. OBSERVATIONS AND CONCLUSIONS 



With, all other parameters held constant, an increase in 
speed drives the design to a very large inside radius while 
reducing the thickness considerably. This trend is so 
strong that a reasonable lower bound on the thickness must 
be given to prevent numerical ill-conditioning. The intro- 
duction of eccentricity increases the effect of rotation on 
the stresses and the optimization produces slightly larger 
values for the design variables. This effect is reduced, 
however, by the constraint on the deflection. 

As the internal pressure is increased, the tendency is 
to produce an optimized radius that is smaller, and a larger 
thickness. This may be explained by the direct effect of 
the internal pressure on the hoop stress which in turn would 
require a larger thickness to support the higher stress. 

The effect of axial load is reflected on the deflection. 
Again, if very small deflections are allowed, the axial load 
effects are reduced. The major parameter that is affected 
is the axial stress. 

A notable characteristic shown by the two variable 
function space diagrams is the tendency of the strength 
constraints to curl back as the radius is increased. This 
signifies the existence of an upper limit on size for a given 
design load. The increase in mass, especially at high speeds, 
increases the stresses considerably leading to violation of 
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the strength constraints. This phenomenon is primarily a 
result of the eccentricity of the mass but also is due to 
the static weight of the shaft itself. 

It is also noted that these strength constraints generally 
run close to each other in an almost parallel manner. The 
dynamic or fatigue constraint maintains a position such that 
coming from the feasible range towards the infeasible area, 
it is activated before the static constraints are encountered. 
In effect, the latter constraints are redundant in rotating 
shafts. However, these are kept in the analysis to allow 
solution of static problems. 

The versatility of the program should not be lost in the 
emphasis on drive shaft applications. With the appropriate 
loadings, one may use it to solve structural problems like 
columns, beams and pressure vessel designs. The objective 
function need not be confined to the volume. For example, 
given a specific initial dimensions, one may minimize the 
Young's modulus to determine whether or not a certain 
material may be utilized; or to determine how low the modulus 
could get before the shaft fails. If the outer dimension 
is critical to the design, one may set a specific value on 
the diameter and let all other parameters vary during 
optimization. 

In summary, numerical optimization provides an efficient 
and effective way of improving shaft designs. The analysis 
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used in this exercise is a first step when compared to the 
boundless refinements in techniques that are possible. 
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V. recommendations for future investigations 



The study has shown the feasibility of using numerical 
optimization techniques in the design of drive shafts within 
the limitations imposed in the analysis. Further studies on 
the same design field may be pursued by eliminating some of 
these limitations. For example, the unfeasible region defined 
by the speed constraint extends only to some discrete distance 
and a supercritical feasible design area does exist. It will 
be interesting to develop optimization techniques to accommo- 
date this disjoint design field and to investigate the pecu- 
liarities this will introduce to the analysis. 

Another investigation may be concerned with shafts with 
variable thickness and radius along its length. The inclusion 
of radial stresses in the calculations, the use of finite 
element methods in three dimensional stress analysis of the 
problem, and vibration considerations are but some of the 
aspects of research that may be useful in future analyses of 
composite materials. 

The challenge of including the unique material character- 
istics of composites in the design will invariably stimulate 
a lot of young minds. However, the handling of these 
direction oriented properties in stress analysis have already 
been investigated [6 through 12], and the tools required to 
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carry out the program are available. COMAND [13] offers a 
simplified composite analysis procedure for panels. Such 
procedures, adapted with the program used in this investi- 
gation could result in a simplified approach to composite 
drive shaft analysis. In this procedure, failure of the 
shaft is assumed if a single ply fails and the entire 
analysis is performed in a ply by ply basis. Since the 
stress elements have to be in a differential level, the 
composite material is assumed to be homogenous and uni- 
directionally isotropic at the differential range. It will 
also be necessary to assume that stresses throughout the 
thickness of the ply is uniform and equal to the stresses 
at the outer element of the ply. Analysis may then be 
performed in the classical procedure, solving for the 
stresses along the major shaft axes. The material properties 
to be used in the calculation must be equivalent values 
derived from the directional properties using transformation 
of axes equations. References 10 and 13 present procedures 
to calculate the corresponding material properties, stresses 
and strains along the major axes of the shaft and the fiber 
orientation coordinates. 

More involved investigations could deal with failure 
modes of the whole shaft instead of a single ply, three 
dimensional analysis of plys including transverse shear 
stresses, delamination problems, thermal effects, and others. 
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The field of composites also opens up another group of 
considerations like the effects of stacking different ply 
orientations at different sequences, unsymmetrical laminates, 
unequal thicknesses of plys and other factors relevant only 
to composites. 
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APPENDIX A 



The FORTRAN program used in this investigation is 
described here. Subroutine ANALIZ provides the basic 
analysis used in the optimization. It reads the initial 
design description and calculates the values of the objec- 
tive function, constraints, and all other parameters 
necessary to solve the analysis problem. COPES-CONMIN 
updates the design to minimize/maximize the objective 
function, iterating until no' further improvement in the 
objective function is possible without violating one of 
the constraints. This process is shown in the flow diagrams. 

The global catalog lists the location, FORTRAN name, 
mathematical symbol and description of the parameters used 
in the optimization. These parameters are contained in the 
labeled COMMOM block, GLOBCM. 

There are at least five input cards to describe the 
initial design. The first card contains the initial 
dimensions — inside radius, thickness and shaft length. The 
second card describes the material used. Young's modulus, 
shear modulus, Poisson's ratio, weight per unit volume, 
yield strength, ultimate strength and compressive strength 
are read here. The endurance limit is internally calculated 
using the endurance limit factors read from the third card. 
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The size factor, CB, is internally updated as the shaft 
diameter is varied during optimization. 

The fourth and fifth cards specify the design and load 
conditions. The parameter, LC, allows the user to choose 
the end conditions applicable to his requirements. Pinned- 
pinned and clamped-clamped analyses are available while 
clamped-pinned and cantilever end conditions are yet to be 
developed. The parameter LCC allows multiple loading design. 
One can thereby design a shaft that is capable of supporting 
LCC different load combinations. 

The magnitude of mass imbalance may be read as a specified 
distance from the center of the shaft, EO, or as a specified 
fraction, EC, of the shaft radius. The safety factor 
generally used in the analysis is also specified in the fourth 
card. 

The fifth card specifies the horsepower, shaft speed in 
RPM, internal pressure and axial load for each loading. 
Additional cards are used up to LCC number of loadings. The 
torque is calculated internally in inch-pounds. The follow- 
ing table summarizes the required information for subroutine 
ANALIZ. 



62 



TABLE IV 



CARD 


INPUT CARDS 
INFORMATION 


FORMAT 


1 


RI, TH, SC 


8F10.2 


2 


E, G, PR, SPWT, I/S 






SU, YC 


8E10.2 


3 


CA, CE, CC, CD, CE, 






CF 


8F10.2 


4 


LC, LCC, EO, EC, SF 


215, 3] 


5 


HP, RPM, P, F 


8F10.2 



Subroutine BEND calculates the bending moment, deflection 
and shear force at the critical cross-sections of the beam. 
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SUBROUTINE ANALIZ (ICALC) 
FLOW DIAGRAM 



ICALC = 1 

READ INITIAL DIMENSIONS 
READ MATERIAL PROPERTIES 
READ DESIGN CONDITIONS 
PRINT ALL INPUT 

ICALC = 2 

CALCULATE GEOMETRIC PARAMETERS 

CALCULATE OBJECTIVE FUNCTION 

CALL SUBROUTINE BEND 

. CALCULATE MAXIMUM DEFLECTION 

. CALCULATE BENDING MOMENTS 

. CALCULATE SHEAR FORCE 

CALCULATE STRESSES AT CRITICAL STRESS 
ELEMENTS 

CALCULATE CONSTRAINTS 
ICALC = 3 

PRINT ANALYSIS RESULTS 



64 



p 




OPTIMIZATION FLOW DIAGRAM 
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TABLE V 

GLOBAL CATALOG OF PARAiVlETERS 



Global 

Location 


FORTRAN 

Name 


Math 

Symbol 


Definition 


1 


R1 


r . 
1 


Inside Radius 


2 


TH 


t 


Thickness 


3 


SL 


a 


Shaft Length 


4 


VOL 


V 


Volume 


5-14 


TQ(IO) 


T 


Torque 


15-24 


P(10) 


p 


Internal Pressure 


25-34 


F (10) 


F 


Axial Load 


35-44 


RDS (10) 


0 ) 


Shaft Speed 


45 


E 


E 


Young ' s Modulus 


46 


G 


G 


Shear Modulus 


47 


PR 


V 


Poisson's Ratio 


48 


YS 


Sy 


Yield Strength 
(Tension) 


49 


YC 


S 

c 


Yield Strength 
(Compression) 


50 


SU 


S 

u 


Ultimate Strength 


51 


RO 


r 

o 


Outside Radius 


52 


SPWT 


y 


Weight Per Unit 
Volume 


53 


AI 


I 


Moment of Inertia 


54 


AJ 


j 


Polar Moment of 
Inertia 


55 


HS 




Hoop Stress 


56 


AST 


^kA 


Axial Stress at A 


57 


ASTN 


(min) 


Minimum Axial Stress 


58-67 


ASC(IO) 


'^xB 


Axial Stress at B 


68-77 


SALT (10) 


0 o ' 

a, a 


Alternating Axial 
Stress 


78 


SMEAN 


% 


Mean Axial Stress 


79 


SST 


Tt 


Torsional Shear 
Stress 
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80 


ssv 


T 

V 


81 


SSMX 


T 

max 


82 


TCRT 


T 

c 


83 


SCSX 


a 

c 


84-93 


VMSI (10) 


a 

m 


94-103 


VMS2 (10) 


o 1 

o 


104-113 


VMS3 (10) 


^mid 


114 


WM 


OJ 

c 


115 


BM 


M 

max 


116 


BMN 


M . 
min 


117 


V 


V 

max 


118 


FC 


F 

c 


119 


SF 


SF 


120-129 


DEL (10) 


6 


130-239 


GC(10,10) 


GC 



Shear Stress due to 
Sheer Force 

Total Shear Stress 

Critical Torque for 
Buckling 

Critical Stress in 
Buckling 

Von Mises (mean) 

Von Mises at x = 0 

Von Mises at Midspan 

Critical Speed 

Maximum Bending Moment 

Minimum Bending Moment 

Shear Force 

Critical Column Buckling 
Load 

Safety Factor 

Deflection 

Constraints 
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APPENDIX B 



The derivation of the equations for the deflection, 
bending moments and shear force at the critical cross- 
sections of the shaft are presented here. Pinnedj ppinned j 
and fixed-fixed end conditions are analyzed. 

1. Pinned-Pinned 




Let y = d sin - 5 - (1.1) 

^ o Z 

2 

dQ = dM (y+e) w 

dQ = pAw^ (y+e) dx (1.2) 

2 

Let = pAco 

dQ = [d^ sin ^ + e] dx (1.3) 
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J 



with the effect of weight: 



dQ = [K^(d^ sin ^ + e) - K 2 ] dx 



(1.4) 



where K 2 = pAg = numerically equal to weight per unit length 



Q 



• r 



- d £ 

dQ = -^+e£] - K2 £ 



(1.5) 



Moment at x = x 



M =-^(£-X) + 
xo 2 o 



I 



(X - x^) dQ 



( 1 . 6 ) 



M = 

X 



= - ? + 


K. 


0 Cl 

^ <TT + l> 

TT 2 


2 


1 


K,e 


2 


- . £. 2 


2 


X 


1 

0 

1 


^2 


2 


K H 


2 


X 


^ 2 ^ 



[§ - " 1 ^ <ir 



TTX 

£ 



(1.7) 



This represents the moment at any value of x. If the axial 
load, F, is included, the total bending moment is 



^TOTAL = ^ ^ 

TT V 

where = Fy = F (d^ sin -^) (1.8) 

Integrating twice: 

Ely = [£ X - 2£x^ + x^] + K, d (— ) ^ sin (1.9) 

^ 24 1 O TT )c 

+ F d (-)^ sin ^ ^ [£^x - 2£x^ + x"^] 

o IT £ 24 *■ 
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At X = l/^, y = 



14 Z 2 ^ K,e2. 

Eld - K, d (-) - Fd (-)^ = — rkir- 

o 1 o IT o 7T 384 



5 £ 

384 



, _ sr 

384 



K^e - K2 



El - - F(^)^ 



] 



( 1 . 10 ) 



This is the deflection when the weight effect is subtractive. 
The maximum deflection occurs when the weight tends to increase 
the deflection. 



max 



51^ 

384 



[ 



K^e + 



K, 



Z 4 Z 2 

El - - F(J)^ 

1 7T TT 



} 



( 1 . 11 ) 



Substituting (1.10) in (1.7) and (1.8), the bending moment 
at X = Z /2 is 



M . 
min 



■'l If* 'So <?>'] 



+ Fd - 






( 1 . 12 ) 



Again the maximum moment occurs when the weight effect is 
additive 



M 



Zs2- 



max 



= [f- + d (-) ] + Fd + 

1 8 max IT ' ■* 






(1.13) 



From Equation (1.5), the maximum shear force is equal 

to I or 



d Z K,e£ K„il 

1 max 1 2 

'^MAX IT 2 2 



(1.14) 
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2 . C 1 amped -Cl amped 

A similar derivation procedure for clamped-clamped ends 
using a shape function: 
d 



O r^ 2tTX, 

y = -^ [1 - COS -^] 



( 2 . 1 ) 



and 



d d 

dQ = [-^ - cos — + e] dx 

d 

Q = ^ dQ = [-^ + e] 



( 2 . 2 ) 

(2.3) 



Moment at x = x 



r 



M = M - ^ X + 1 (x - x) dQ 

0 2 J O 

O 



K, 2. d 

_ A/t 1/0 

“ ^o ■ 2 ^2 



+ e) + r° 



X d X d 

r O o o o 27TX 

[ = ~ ^ COS Q 



xd xd^ 

+ X e - — cos — ^ - xe] dx 
O 2 ^ ^ 



(2.4) 



M at any value of x: 

K il d K, d ^ 

Ely" = M = M ^ (^ + e) X + -^ (^ + e) x 

X O 2 2 2 2 



Kid H 

1 o /T 2 ttx^ 

„ 2 ' —> 

8tt 



(2.5) 
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Integrating once; 



d 32. 

Ely' = X - ^ + e) (^ - 2ji) 



K, d £ X 
1 o 

2 

8tt 



K,d V 
1 o 

167T^ 



Sin IS . c, 



y' = 0 at X = 0, 1/2 and = 0 

«o = "^1 ^ A 

Sir 

Integrating once more; 



r-T d 4 3„ 2„2 

^ = f_o + e) (- ^ ^ ^ ) 

^ 2 '24 12 24 ^ 



d a 

O 

32tt' 



2 ttx , ^ 
cos — + C2 



y = 



0 at X = 0, ^2 ~ 



d 1 
o 

32 t \‘ 

2 






(1 - 



cos ^4^) (2.6) 



y = 



d at X = Jl/-,; 
o 2 



|iao= ("s.e, 4 , 4 . 



d I 
o 

167t^ 



Eld 
c 

Kn 



4 4 

d^il^ ^.4 d^Jl^ 

o e t o 



768 



384 



167T 
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(2.7) 



d 



o = 



K^eJl 



C 4 

384EI + 24) 

TT ^ 



With axial load F; 



Fd 



Moment due to F = M^. = M^ - F^ = M^ - (1 - cos — — ) 



2ttX\ 

£ 



Integrating once: 



( 2 . 8 ) 



' 0 9 

Elyp = X - ^ (X - sin + Cj_ 



0 at X = 0, C, = 0 



. Fd 

yp =0 at X = e, = — 



where M = Moment at the ends due to the axial load 
o 

Fd 2 Fd 2 „2 ^ 

- O ,x . O rX J- 2. 2tTXt , _ 

2 ^ 2 ^ 2 ^2 a 2 ^ ^2 

4it 



Yp = 0 at X = 0, = 



Fdo£ 

8 u 2 



Fd £ 

o /n ^ 2 tTXv 

Yr- = — o — (1 - cos -!r-) 



F 2 

8tt^EI 



£ 



(2.9) 



At X = 1/2 J 



Fd £ 

'f = ~1~ 

4tt El 
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From Equation (2.6) 



K, d 2 « 2 - 

= Ef <f4> ® <54* 



d I' 
o 

4 



/ 1 2ttXx t 

(1 - cos -5—) ] 



At X = a/ 2 '• 



4 4 

£ d Kt d £ 

y = r_2. + 0 ] 4- _i 2 

384EI '■2 T . 4^^ 

16it EI 



■^O ' i'p 

= (J!iL + ''1“' 



4u^EI ' 167 t"EI 



4 

£ e 

^ d 4 - i 

4„,'' o 384EI 



q 2 K £^ 24K £^ 

384EI d = ^ 

o J2. 



TT 



d = 
o 



2 TT 

4 

K^e£ 



^ ) d,+ el 



384EI - Kj^ l^ ■ 96^ 

IT TT 



When the effect of weight is considered: 



^o = 



4 4 

K^e£ ± K2£ 

384EI - 0.74638K^£^ - 96(|)^F 



where K 2 = pAg, numerically equal to weight per unit 



(2.10) 



( 2 . 11 ) 



( 2 . 12 ) 



length. 
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